WIS-94/3/Jan-PH 
January 1994 



o\ '■ Phenomenology of models with more than two Higgs doublets 

Yuval Grossman 
en 
(N 

Department of Particles Physics 

The Weizmann Institute of Science 
^^ , Rehovot 76100, ISRAEL 

m 

o 

^ ■ We study the most general Multi-Higgs-Doublet Model (MHDM) with Natural Flavor Con- 

Q-i, servation (NFC). The couplings of a charged scalar Hf^ to up quarks, down quarks and charged 

O ' leptons depend on three new complex parameters, Xi, Yi and Zj, respectively. We prove relations 

ZJ . among these parameters. We carry out a comprehensive analysis of phcnomenological constraints 

1-^ ' on the couplings of the lightest charged scalar: X, Y and Z. We find that the general MHDM may 

differ significantly from its minimal version, the Two-Higgs-Doublet Model (2HDM). 
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Abstract 



1 Introduction 

The Higgs sector of the Standard Model (SM), consisting of a single Higgs doublet, has not yet 
been experimentally tested. The possibility of an extended Higgs sector is definitely consistent with 
experimental data. The simplest extensions are models with several Higgs doublets, implying the 
existence of charged scalars. The simplest of these is the Two-Higgs-Doublet Model (2HDM), which 
has received substantial attention in the literature [Q| . However, very little is found in the literature on 
models with more than two Higgs doublets and, in particular, there exists no comprehensive analysis of 
these models. It is the purpose of this work to study the general Multi-Higgs-Doublet Model (MHDM), 
and perform a comprehensive analysis of the charged Higgs sector of the model. In the lack of direct 
experimental knowledge, we study indirect effects of charged scalars in order to place constraints on 
this sector. In particular, we consider charged scalar contributions to low-energy processes: neutral 
meson mixing, rare K and B decays, Z boson interactions and CP violating effects. The results from 
this combined analysis are used to find the region in parameter space which is consistent with all data. 
We then discuss the differences between the MHDM and the 2HDM. Throughout this work we assume 
that there are no significant cancellations with other types of new physics. 

Our work is organized as follows. In chapter 2 we describe the MHDM and prove general relations 
between the model parameters. In chapter 3 we discuss the relevant processes, and list the formulae 
and numerical data used in our analysis. In chapter 4 we combine all the constraints and find the 
allowed range of the model parameters. In chapter 5 we discuss the differences between the MHDM 
and the 2HDM. Finally, chapter 6 contains a summary of our work. 

2 The Model 

2.1 Introduction 

There are basically two major constraints on any extension of the Higgs sector of the SM. First, it 
is an experimental fact that p = m'y^r / {rn^ cos"^ 6\y) is very close to 1 [p. It is known [|l| that in a 
model with only Higgs doublets, the tree-level value of p = 1 is automatic without adjustment of 
any parameters in the model. Second, there are strong experimental limits on flavor-changing neutral 
currents (FCNC's). In the SM, tree-level FCNC's are absent because fermion mass matrices and 
Higgs-fermion couplings are simultaneously diagonalized. In general, this ceases to be true in a model 
with a non-minimal Higgs sector. An elegant way to avoid this problem is based on a theorem by 
Glashow and Weinberg |y] called Natural Flavor Conservation (NFC): tree- level FCNC's mediated 
by Higgs bosons will be absent if all fermions of a given electric charge couple to no more than one 
Higgs doublet. This can be achieved by imposing extra symmetries. If we adopt this mechanism, the 
Higgs couplings to fermions are constrained but not unique. There are five possibilities to couple the 
Higgs doublets to the known three types of massive fermions (up- type quarks, down- type quarks and 
charged leptons.). Summary of these possibilities is given in Table ||. For Model III at least three 
Higgs doublets are needed, for Model I one is enough, while for the other models two are sufficient. 
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d (down-type quarks 
u (up-type quarks) 
e (charged leptons) 

Table 1: Summary of all the possibilities of MHDM. The numbers in the table show which Higgs 
doublet couples to which fermion type. 

It does, of course, make a difference if the number of Higgs doublets is larger than the minimal one. 
For a general MHDM the Yukawa interactions are given by 

-Cy = QL.^F,pdDR^ + Ql^F^j^uUr^ + LLApeER^ + h.c. (1) 

Where left-handed quark doublets are denoted by Q^. , and left-handed lepton doublets by L^. . Right- 
handed quark singlets are denoted by Dr^ and Ur. , and right handed charged lepton singlets by Er. . 
The sub-index z is a generation index (i = 1,2,3). Higgs doublets are denoted by $j (where j runs 
from 1 to n, the number of Higgs doublets), and $j = i(T2<l>*. Sub-indices d, u and e denote the Higgs 
doublet that couples to the down-type quarks, up-type quarks and charged leptons, respectively. F^ 
and F are general 3x3 Yukawa matrices and we can choose a basis where one of them is real and 
diagonal. Since we consider massless neutrinos, F can be chosen real and diagonal. In a general 
MHDM, with n Higgs doublets, there are 2n charged and 2n neutral scalar fields. After spontaneous 
symmetry breaking two charged fields and one neutral field become the would-be Goldstone boson 
"eaten" by the W and the Z in order to acquire their masses, and only 2(n — 1) charged and 2n — 1 
neutral physical scalars are left. The Yukawa interaction of the physical charged scalars with fermion 
mass eigenstates is given by 0] 

n 

Cp = (2V2Gf )i/2 Y^(x,UlVMdDr + YdJRMuVDL + Z,NlMeEr)H^ + h.c. (2) 

i=2 

Mjj, Mu and Me are the diagonal mass matrices of down-type quarks, up- type quarks and charged 
leptons, respectively. H^ denote the positively charged physical scalars. Left-handed neutrino fields 
are denoted by A'^^, and the CKM matrix by V. Xi, Yi and Zi are complex coupling constants 
that arise from the mixing matrix for charged scalars. An important feature of the charged scalar - 
fermion interactions is that they are proportional to the fermion masses. Hence, effects are stronger 
in processes involving heavy fermions, real or virtual. 



2.2 Relations between Xi, Yi and Z, 

Xi, Yi and Zj arise from the mixing matrix U, which rotates the charged scalar interaction eigenstates 
into the mass eigenstates y] 
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where G^ is the would-be Goldstone boson and H^^ are the physical charged scalars. U is an n x n 

unitary matrix that can be parameterized as the generalized CKM matrix with n{n — l)/2 angles and 

(n — l){n — 2)/2 phases Q. Since G"*" is the Goldstone boson that gives the W~^ its mass in the same 

way as in the SM, we obtain 

1 



G-' = -Y.-^'^t, (4) 

where Vi =< ^i > and v = 246 GeV is the SM vacuum expectation value (vev). Using eq. (0) we get 



Vi = UiiV, 
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Using these definitions and the unitarity of U, we obtain the following relations: 



Relation 1 
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Y^YZ* = l (for CD, / CD,). 
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Relation 3 Xi,Yj and Z^ {i,j,k S 2..n) cannot be simultaneously large (more than \/n — 1), unless 
there is at least one Higgs doublet that does not couple to fermions. 



These relations are proved in appendix A. 

Special cases of our results are known for 2HDM and 3HDM. In 2HDM there is only one charged 
scalar, so we drop the sub-indices of Xi, Yi and Zj. For Model II [|l| 

Z = X = tan/3, y = cot/3, tan/3 =—, (9) 

and the three relations are clearly fulfilled. For 3HDM P] 

X,Y^ = l-X,Y^, Yl + Y:i = '^A±^^ Xl + xl = 'S^, (10) 

which are examples of relations |l| and |^. 
A few points are in order: 

1. An important relation between X and Y in 2HDM (Model II), namely XY* = 1, is not retained 
in the framework of a general MHDM. Consequently, processes involving this combination can 
be enhanced or suppressed in MHDM compared to 2HDM. A similar statement holds for the 
combinations XZ* and YZ*. 

2. The most general features of MHDM can be realized in 4HDM. In 3HDM {Model III) the 
couplings cannot be simultaneously large. This constraint is removed in 4HDM. However, in 
most of the processes that we study, not all couplings are involved, and 3HDM is general enough. 

3. From relation |l| we see that if all charged scalars were degenerate, the model practically reduces 
to 2HDM, namely, the relations XY* = 1, XZ* = — 1 and YZ* = 1 hold. Consequently, in order 
to exploit the most general MHDM, we assume that one of the charged scalars is much lighter 
than the others, and that all the heavy charged scalars effectively decouple from the fermions. 
Then, we consider in the calculations the single light charged scalar, and drop the sub-indices 
ol Hi, Xi, Yi and Zj. 

2.3 CP violation 

Within the framework of MHDM, CP violation could arise in charged scalar exchange if there are at 
least three Higgs doublets [g]. Then, there are two mechanisms which give rise to CP violation (for a 
review see [^): 

1. Spontaneous CP violation (SCPV). CP symmetry is enforced at the Lagrangian level, but broken 
through complex vev's. The requirement of SCPV forces <5ckm = [^. In this case, CP 
noninvariance arises solely from charged scalar exchange. However, experimental data exclude 
this possibility |^, ^. 

2. Explicit CP violation. CP is broken by complex Yukawa couplings and possibly by complex 
vev'sfj. In such a case, CP violation can arise from both charged scalar exchange and W^ 
exchange. 



^ Actually, CP can also be broken through complex scalar couplings. We will not discuss this possibility. 



Figure 1: The Feynnian rules for the charged scalar - fermions interaction within MHDM. 

In both cases, CP violation in the Higgs sector is manifest in phases that appear in the combinations 
XY* , XZ* and Y Z* . We emphasize that the requirement that the charged scalars are nondegenerate 
is mandatory for CP violation from the Higgs sector. If the charged scalars were degenerate, XY* is 
replaced by Y^^=2^i^i' = 1) which is real. 

2.4 Summary 

The model that we study is the SM with extended Higgs sector. Its Higgs sector contains more than 
two Higgs doublets with NFC. We focus on the charged scalars, and assume that all but the lightest 
of the charged scalars effectively decouple from fermions. The couplings of the lightest charged scalar 
to fermions are given by: 

£+ = {2V2GFf''^{XU lV MdDr + YUrMuVDl + ZNlMeEr)H+ + h.c. (11) 

This interaction yields the Feynman rules given in fig||. In this model there are several new parameters: 
The mass of the lightest charged scalar: mn, and the coupling constants: X, Y and Z. These couplings 
are arbitrary complex numbers. It is the subject of the coming chapters to find the experimentally 
allowed ranges for these parameters. 

3 Phenomenology 

In this chapter we discuss various processes where a virtual charged scalar may contribute significantly. 
From the experimental results for these processes, we are able to constrain the parameters of the model. 
Wherever possible we use the experimental data at the 95% CL. For theoretical uncertainties we use 
estimates and cannot assign a confidence level to the errors. Before discussing the processes, we list 
values of quark sector parameters that we use, and give constraints on the charged scalar parameters 
that are obtained from direct measurements and perturbativity. 

3.1 Values of quark parameters 

For quark masses we use pOt| 

m„ = 5MeV, m^ = 9 MeV, m^ = 180 MeV, m^ = 1.4-1.5 GeV, mfe = 4.6-5 GeV. (12) 

For the top quark the situation is different. First, it is the only quark that has not been observed. 
Thus, there is large uncertainty in its mass. Second, since it is the heaviest quark, its coupling to the 



charged scalar is the largest. The direct search for the top at CDF gives mt > 113 GeV [11|, while 



the electroweak precision measurements at LEP give mt < 185 GeV [12|. Within the framework of 



MHDM both bounds are relaxed. If mn < nit, the top may escape detection at CDF [^]. Additional 



loops that include the extra scalars may affect the electroweak precision measurements [14|. Then, 



the bounds on rut are related to the MHDM parameters. Thus, whenever we use some values of the 
MHDM parameters we use only values of rrit that are consistent with them. 
For the values of the CKM matrix elements we use |15, ^, 17] 



\Vud\ = \Vcs\ = \Vtb\ ^ 1, (13) 

\Vcd\ = \Vus\ = 0.22, 

\Vts\ = \Vcb\ = 0.040 ±0.007, 

\Vub/Vcb\ = 0.08 ±0.03. 



|Fid| =0.004-0.016. (14) 



Unitarity constraints on \Vtd\ give 



The allowed ranges for \Vtd\ and 5ckm are different in the SM and in MHDM ||TM|. Furthermore, they 
are related and depend on rrit- 

For the masses of the physical particles we use the values quoted by the Particles Data Group 
(PDG) @. 

3.2 Constraints on scalar parameters 

For the mass of the lightest charged scalar we use |jl^ 

ruH > 43 GeV. (15) 

There is no upper bound on niH'- The SM Higgs boson is bounded by unitarity to be lighter than 
about 1 TeV, but there is no similar unitarity bound on the charged scalar of MHDM (for a discussion 
of this point see |]l|). 

We can get perturbativity bounds on \X\, \Y\ and \Z\: We restrict ourselves to the range where 
perturbation expansion is valid. Following [EO] we find 



rrit 

nib 
|Z| < 55^^ » 340. 

There is no lower bound on \X\, \Y\ and \Z\. When all of them vanish, the model practically reduces 
to the SM. 

3.3 B — B mixing 

The mass difference in the neutral B meson system (Ams) has been measured and it fits the SM 
predictions |21, ^, ^. However, there are large uncertainties in the calculation (mainly in the values 



Figure 2: The diagrams that contribute to B — B mixing, (a): The SM diagram, (b): The two-Higgs 
diagram, (c): The mixed diagrams. The diagrams where the intermediate bosons and fermions are 
exchanged are not shown. 

of mt, Vtd and fs) that leave a lot of room for non-SM contributions. The MHDM extra contributions 
to the mixing come from box diagrams with intermediate charged scalars (see fig. |^. 

In the calculation we neglect external momenta but keep all quark masses. We use the Vacuum 
Insertion Approximation (VIA) which is a good approximation in this case ||2^. For intermediate 
quark only the top quark is important 

MuiBg) = -^ml,r]BiVt*gVtbf[Iww + Ihh + Irw] q = d,s. (17) 

The SM contribution (fig. 0(a)) is /vi/H/[p^], while box diagrams with two charged scalar propagators 
(fig. §(b)) or one charged scalar and one W propagator (fig. |2|(c)) give Ij^jj and Ijj^y, respectively. 
r]B ~ 0.55 is a QCD correction factor pi]. The QCD corrections to the MHDM diagrams have been 
calculated in the leading logarithmic approximation [^], and only for the terms with the same Dirac 
structure as the SM diagrams. However, following |20| we apply the same QCD correction factor both 



to the SM and the MHDM contributions. In what follows we present only the most important terms: 

Iww = xt4.Io{xt)VLL, (18) 

Ihh = xt[yt\Y\^Ii{yt)VLL + ytyb{XY*)Hh{yt)SLL + ybyq\xfh{yt)VLLi 

IhW = Xt[xt\Y\^{2h{xt,XH) -8h{xt,XH))VLL + '2xtXb{YX*)h{xt,XH)SLL], 

where Xg = m'^/m'yy and yg = m^/m^j. The matrix elements, Vll and Sll-, and the loop integrals, 
/j, are given in appendix B. Setting rrn, and m^ to zero we get the known results p^, p^, while for 
the case where only m^ is set to zero our results agree with those of Q. We use Am-s = 2|Mi2| 
and g |2|, 



Xrf = Ami3/rB = 0.67 ±0.11, /b = 190 ± 50 MeV, rg = 1.49 ± 0.03 ps. (19) 

We conclude that large contribution to B — B mixing may arise for large Y or large XY* . For 
the B system the contribution from the term proportional to |X| is small, even for X as large as the 
perturbativity bound. However, for the Bg system its contribution can be large. 

3.4 K — K mixing 

There are two well measured quantities related to K — K mixing: the mass difference /\mK and the 
CP violating observable e (for a review see |l^ ) . Taking into account the bounds from B — B mixing, 
we find that charged scalar mediated diagrams may give only a small contribution to Atjik ||1|- The 
situation is different for e |22| 

\e\ = '4^. (20) 



e is a very accurately measured quantity, but there are large uncertainties in the SM calculation 
(mainly from mt, Vtd and Bk)- In MHDM there are two new sources that can contribute to e: The 
imaginary part of the charged scalar box diagrams, and enhanced long-distance contributions. The 
short-distance diagrams are similar to those oi B — B mixing (fig. ^ with the replacement b ^ s. 
For the SM box diagrams we consider both charm and top intermediate quarks, while for the charged 
scalar box diagrams only the top quark is important. We keep only the terms that are proportional 
to the internal quark masses: 

\esD\ = ^J^^2^w^KflBK Im[r + /* + 1% (21) 

with 

r = r]ixXh{xc). (22) 

/* = ri2Xt\][lQ{xt) + yt\Y\^h{yt)+Xt\Y\^{2h{xuXH)-^h{xuXH))l 
I"*- = 2r]3XcXtXch{xc,xt), 



where A^ = VqgV*^, and r]i ~ 0.85, ri2 ~ 0.62 and r/3 ~ 0.36 are QCD corrections factors [25|. Due to 
the lightness of the valance quarks, the terms that depend on their masses are neglected. This implies 
that the term that may have introduced CP violation from the Higgs sector, namely, the term that 



depends on XY*, is very small and cannot contribute significantly to e. Our result agrees with |29|. 
For the numerical values we use P, E2| 



|e| = 2.26 X 10"^ fx = 161 MeV, Bk = 2/3 ± 1/3. (23) 

The long-distance contributions to e are negligible in the SM but in MHDM, where they depend 



on the phases of the Higgs sector, they may be important (see [22, ^, ^). The largest contribution 

,2 im(xy*) 



comes from the intermediate rj' |31, 22 1 



I ELD I -0.02 GeV 



rrir 



In ^-1.5 
mj. 



(24) 



There is an additional CP violating parameter, e'/e. There is no unambiguous evidence for e'/e 7^ 
0]. Furthermore, there are large uncertainties in its calculation. Therefore, at present, it is not 



possible to extract useful constraints from it. A discussion can be found in [|y, 18]. 

We conclude that large contribution to e may arise for large Y through the short-distance contri- 
butions, or for large Im(Xy*) through the long-distance contributions. 

3.5 D — D mixing 

The mixing in the D system is different from other neutral meson systems: 

1. The intermediate quarks in the box diagrams are down- type. Consequently, the box diagrams 
contribution to the mixing is smaller than in mesons where the intermediate quarks are up- 
type Hi. 



2. D — D mixing has not been observed, and only an upper bound has been estabhshed |^, 

AruD < 1.5 X 10"^^ GeV 95% CL. (25) 

3. The long distance contribution to the mixing is probably dominant over the SM box diagram. 



It is estimated to be about two orders of magnitude below the experimental bound [33| or even 
smaller [34|. 



Because of the above reasons, D — D mixing constrains neither the SM parameters |15|, nor (as we 
find) the MHDM parameters. 

The diagrams that contribute are similar to those oi B — B mixing (fig. |2|) with the replacements 
b ^ c, d ^ u and t ^ s,b. In the calculation we neglect the SM contribution and the mixed diagrams, 
while from the two-Higgs box diagram we take only the term proportional to |X| 

G ^ 

^12 = -^^rnl^fD^DlHH^ (26) 

with 

Irh = \M'^XH[>?sylh{vs) + >^lylh{yb) + '^hKyhyshivh^ys)], (27) 

where \q = VcqV*g. We can simplify the expression using 

hiys) = hiyb) = hiyb,ys) = h (28) 

and get 

Mi2 = ^m2^/|>mz)|X|^XH[A2y2 + x^y^ + 2XbXsybys]- (29) 

Setting Aft to zero, our result agrees with (2^, while setting nis to zero, our result agrees with pO|| . 
For the numerical values we use IQ, ^ 

AniD < 1.5 X 10~^^ GeV, fo = 170 ± 30 MeV. (30) 

We conclude that large contribution to D — D mixing may arise for large X. 

3.6 Z ^bb 

The decay Z ^ bb has been measured at LEP. Within MHDM, extra contributions arise via virtual 
scalars in the vacuum polarization diagrams and the vertex corrections. Recently, a new method 
of analyzing the electroweak precision data was proposed [^5|. In this method the non-SM effects 
manifest themselves via four observables: ei, €2, £3 and e^. We are interested in et, which is closely 
related to the vertex correction to Z — > bb. By comparing the experimentally measured quantity, 
et,, to the theoretical calculation of the vertex correction, V;,, one can put bounds on the MHDM 
parameters [|^ 

eb = 2.3V6. (31) 



Figure 3: The additional MHDM diagrams that contribute to b 
replacing the charged scalar propagators with a VF-boson ones. 



S7. The SM diagrams obtained by 



Vfc is defined in |37]: 



V, 



V5Af + vH 



b ' 



with 



Vpf 



20a 



mt 



rrir 



(32) 



(33) 



V 



H 



C ■ FL{mt)\YYml 



where C is given in appendix B and Fl can be found in |^7[. For the numerical data we use |^5|] 

€6 > -1 X 10-2 95% CL. (34) 

There are also terms in V^ proportional to Im^Xp. Some of these terms depend on the neutral 
Higgs parameters, but appear with a sign opposite to the terms that depend on the charged Higgs 
parameters [^]. Thus, cancellation may occur and a bound on X cannot be obtained. This also 
applies to the decay Z -^ rf. There, the MHDM contribution is proportional to \mrZ\^, but the 
terms that depend on the neutral Higgs parameters may cancel the terms that depend on the charged 
Higgs parameters, and a bound on Z cannot be obtained. 

We conclude that large contribution to Z — > 66 may arise for large Y or large X. However, a bound 
can be obtained only on Y . 

3.7 6-^s7 

Recently, the CLEO collaboration has given an upper bound on the radiative B decay [39| 

BR(6 -^ 57) < 5.4 X 10"^ 95% CL, (35) 

which is only about a factor of 2 above the SM prediction. Uncertainties in the calculation come 
mainly from the value of mt |4C]. The diagrams that contribute to 6 ^ 57 are given in fig |3[ The SM 
calculation was done in [41| and the MHDM calculation in [^, ^, ^, 44 1: 



BR(6 ^ 57) = C\m + Gwixt) + {\Y\^/3)Gwiyt) + iXY*)GHiyt 



where 



C = 



3arif BR{B -^ XJu) 
27rFp,(m2/m2) 



3 X 10" 



(36) 



(37) 



Fpsr^ 



0.5 is a phase space factor given in appendix B, 771 ~ 0.66 and 772 ~ 0.57 are QCD correction 



factors [44|. The expressions for the G-functions are given in appendix B. The SM result can be 
obtained by setting X and Y to zero in (pq). In the full calculation [^] a term proportional to |X| 
appears, but it is suppressed by (ms/mi,)'^ and we can safely drop it. 

We conclude that large contribution may arise for large Y or large XY*. 
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Figure 4: The diagrams that contribute to K'^ -^ n'^vu. (a) The Z penguin diagrams. The sohd square 
represent the effective sdZ vertex which induced by diagrams similar to these of 6 — > ,37 (fig. 0). (b) 
The SM box diagram, (c) The two charged scalars box diagram, (d) The mixed diagram. The second 
mixed diagram is not shown. 



3.8 K+ 



TT ' uiy 



The branching ratio of the decay K^ -^ tt^vu has not yet been measured. The current bound is [46| 



BR{K+ -^ TT+VT?) < 6.8 X 10"^ 95% CL, (38) 

is about two orders of magnitude above the SM prediction. The calculation suffers from uncertainties 
in rrit and Vtd (for a review see ^^). The diagrams that contribute to this process are given in fig. |^. 
In the SM, the penguin and the box diagrams are of the same order, and intermediate charm and 
top quarks are important. In diagrams that involve charged scalars the intermediate top quark is 
dominant. Since the external quark masses are small, we consider only terms proportional to internal 
quark masses. 

The SM diagrams were calculated in [41|. The MHDM result is obtained from the 2HDM calcula- 



tion [pO| , |l8| , ^ by multiplying the two charged scalars box diagram by \YZ\'^ and the mixed diagrams 
by Re{YZ*) 



\Vu 



with 



C ^ ^I^BR(i^+ ^ vrOe+z.) ^ 6.9 x 10"^. (40) 

The SM contribution is / , the two charged scalars box diagram is / , the mixed diagrams are 
/ and the charged scalar penguin diagrams are / 

I^^ = r]AcDsM{xc,Xi) + XtDsM{xt,Xi), (41) 

jHH ^ \YZ\^XtxtyfMyt,yi), 

jWH ^ Re{YZ*)XtXtXi{l7{xt,Xi,XH) + l8ixt,Xi,XH) ), 

I^H ^ \Yf\tDzH{xuyt), 

where Xg = VgdV*^ and the D-functions are given in appendix B. rjc ~ 0.7 is a QCD correction 
factor ^^, while QCD corrections to the diagrams with an intermediate top quark are small |p9|| . 
While for the penguin diagrams the MHDM contribution is independent of the final lepton, the 
charged scalar box diagrams depend on the internal lepton mass and it may be important only for 
internal r. In 2HDM (Model II) YZ* = 1 and the charged scalar box diagrams are negligible. Our 
results disagree with those of ^^l where it is claimed that the charged scalar box diagrams are not 
important even within MHDM. We believe that the reason for this is that they assume that all the 
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Figure 5: The diagrams that contribute to B ^ Xtu^-. (a) The SM diagram, (b) The charged scalar 
mediated diagram. 

coefficients in the Higgs potential are of the same order. Then, \YZ\ = 0(1) and the effect is indeed 
small. 

A similar calculation applies to the decay B -^ Xguv. However, due to the lack of experimental 
data, no constraints on the parameters can be obtained. 

We conclude that large contribution to K^ -^ n^vu may arise for large Y or large Y Z* . 

3.9 B -> XTPr 



The branching ratio for B — > XtUt has been recently measured by the ALEPH collaboration [50| 
which leads to the upper bound 

BR(S -^ Xrvr) < 4.00% 95% CL. (42) 

The uncertainties in the calculation are mainly from m^ and rric- The diagrams that contribute to 
the decay are given in fig. |5l In the calculation we neglect charmless final states, since they are 
highly suppressed by CKM elements. A recent SM calculation was done in [£^]. The MHDM result is 
obtained from the 2HDM calculation |52, ^, 54] by replacing |X| with XZ* 



BR{B -^ Xrur) = BR^^\B -^ XTUr){l + j\R\^ - D ■ Re(i?)), (43) 



where the SM result is |51] 



BR^'''{B -^ Xrvr) = (2.30 ± 0.25)%, (44) 

and 

^ _ mrm^XZ* ^ _ ^ rrir F^f(mc/mb,m^/mfc) 

ni^j ' mb Fps{mc/mb,mr/mb) 

The phase space functions Fps and F*"* are given in appendix B. For the ranges of masses that we use 
D ~ 0.43 lb 0.01. We disagree with [p^ j about the numerical value of D. Due to this disagreement the 
bound that we obtain is weaker than the bound obtained by [^]. In 2IIDM (Model II) the interference 
term always reduces the rate, while in MHDM it may also enhance it, depending on the relative phase 
between X and Z. Nevertheless, in order to get bounds on \XZ\ one should take the minimal value 
for the branching ratio, which corresponds to aTg{XZ*) = 0. Thus, the bound on \X\ in 2HDM is 
the same as the bound on \XZ\ in MHDM. 

We conclude that large contribution to B ^> Xrv-r may arise for large XZ. 

3.10 e— /i universality in r decays 

e-jJL universality is confirmed in leptonic decays of pions and kaons, and in r decays (for a review 
see p5|). The SM and MHDM diagrams are similar to those of -B — > XtVj- (fig. |5|) with the replace- 
ments: 6 — > r, c — > z^T- and r — > £ (£ = e,/i). At first glance, it seems that a useful bound on the 
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MHDM parameters cannot be obtained since the final leptons are very light. However, in purely lep- 
tonic decays both the theoretical calculation and the experimental measurement are clean. We then 
require that the charged scalar contribution does not exceed the difference between the SM prediction 
and the experimental result. The calculation in the case of muon decay was done in ||5^ and for r 
decay in 



I 00 

BR(r -^ ^J.yrV^l) 
BR(t — > evT-Ve) 



B? B? 

Fps{ml/ml)[l + — -D-R]^ 0.9726[1 + -- - 0.U5R], 



where 



R 



mr'm^\Z\ 



D 



,m^ F;f (ml/ml) 



(46) 



(47) 



Fp" , the phase space function for the interference term, is given in appendix B. We use the recent 



numerical values 



which leads to the upper bound 
BR(r -^ /ii^ri^^) 



< 1.035 95% CL. 



(48) 



BR(r -> eVr^e)Fps{ml/ml) 
We conclude that strong deviation from e-// universality in r decays may arise for large Z. 

3.11 Electric dipole moment of the neutron 

An electric dipole moment (EDM) of an elementary particle is a manifestation of CP violation (for a 
review see [^]). In our discussion we concentrate on the EDM of the neutron (NEDM), Dn- Although 
there are very large uncertainties in the calculation, we know that the contribution to the NEDM from 
the electro-weak sector of the SM is much smaller than the current experimental bound [0] 



DJ < 1.2 X 10"^^ e cm 95% CL. 



(49) 



A large NEDM can be generated within MHDM in many ways (for a review see [30|). In our discussion 
we concentrate on the two that involve charged scalars: the EDM of the down quark, and the three 



gluon operator. The contribution to the NEDM from the EDM of the down quark is given in |59|. We 
neglect the contribution from the internal top quark and using rric <C mn we get 

V2Gf 



Dr, 



97r2 



-rud lmiXY*)\Vcdr (HVc) + 0.75) y^. 



For the three gluon operator we use the naive dimensional analysis result 

4.V2CGFg^, 



Dn 



{AttY 



-lm{XY*)h'{yt), 



(50) 



(51) 



where 



h'{y) 



{i-yY 



"^y-y 



21n(y) 



We disagree with the numerical factor of the interference term in eq. (8) of 
Fps^ /Fps instead of muhiplying by this factor. 



(52) 
It seems to us that they divide by 
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(, the QCD correction factor, depends sensitively on the scale [61] 



( = 10"^ -10"^ 
We conclude that large contribution to the NEDM may arise for large lm.{XY*). 



(53) 



3.12 B -^ i+i- 

The branching ratio for B -^ i'^i~ has not been measured, but there is an upper bound for i = e, fi. 
The diagrams that contribute to B ^ i'^i~ are similar to the diagrams that contribute to K^ -^ n^vv 
(fig. §). The photon penguin diagrams vanish because of electromagnetic current conservation. The 
charged scalar penguin diagrams contribute universally for all leptons. However, the charged scalar 
box diagrams depend on the lepton mass. Thus, it is important only for B -^ t^t^ . We emphasize 
that in 2HDM the box diagram is negligible since it is proportional to \YZ\'^, which is 1 in Model II, 
or 0(1) in Model I. A recent SM calculation was done in |£2| and the 2HDM in ||6^. We calculated 
the box diagram. We neglect terms proportional to external quarks masses, and the QCD correction 
which is small [H^]. We get 



BR{B -^ i+r) = BRsm{B -^ i+r 






1 2 



where 



BRsm{B 



p+f-'] 



TB 



327r5 



4m| 



flmBmj\ll _2 



m 



B 



G'''\xt)Xt 



(54) 



(55) 



Xq = ^qdVqh a^d G^ and G^^ are the charged scalar penguin and box diagram contributions, respec- 
tively: 



Gf 



ytXt\Y\ J{yt), 



(56) 



GL = -xtyei\YZ\^Ig{yt) + 2ReiYZ*)ho{xt,XH)). 



GsM, J and Ij are given in appendix B. Our results agree with |Q, but disagree with ^^ about the 
Dirac structure of the matrix element, and the sign of the first term of /g. However, since they work 
within the framework of 2HDM, where the box diagram is negligible, this disagreement does not have 
numerical significance. 

A few points are in order: 



1. The same calculation applies to Bg decays with {d 



s). Consequently, the branching ratio for 



Bg decay is larger by approximately li^^'^ 



2. It was noted |64, |6^ that large effects may arise due to neutral scalar penguin diagrams, even 
in 2HDM. This contribution depends on neutral Higgs sector parameters and could also be 
small. We note that the contributions from the charged scalar and from the neutral scalars are 
independent. 
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3. The process Ki^ -^ ^^^^ gets contributions from diagrams similar to i? ^ i'^i^ . Since the muon 
is hght, the box diagram is neghgible. In addition, the long distance contribution dominates. 
Thus, it is difficult to estimate the charged scalar penguin diagrams contribution. We do not 



discuss this process. For a detailed discussion in the framework of 2HDM, see |2C]. 



The processes B -^ Xsi~^i gets contributions from diagrams similar to B ^ l^i , with the 
addition of photon penguin diagrams. In 2HDM the charged scalar box diagrams are negligi- 



ble |42, ^. However, in MHDM they may be important for final r leptons. Then an interesting 
situation arises: The Z penguin depends on Y, the photon penguin depends on Y and YX*, 
and the charged scalar box diagrams depend on \YZ\. The process B -^ XsT~^t~ then depends 
on all the charged Higgs parameters. Furthermore, the neutral scalar penguin Q, ^] may also 
contribute. However, due to the poor experimental data, no bounds can be obtained from these 
decays. 

We conclude that large contribution to B ^ i'^i~' may arise for large Y, and for B -^ t^t^ also 
for large YZ. 

3.13 B ^ iu 

The branching ratios of the decays B ^ iu {i = fJ.,T) have not yet been measured, and the current 
90% CL bounds ^ 

BR{B -^ ^JLu) < 2.0 X 10-^ (57) 

BR{B -^ Tv) < 1.2%, 

are about two orders of magnitude above the SM predictions. The calculation suffers from uncertainties 
in fs and Vub- The diagrams that contribute to the decays are similar to those of i? — > XtVt (fig. |5|), 
with the replacement c ^ u. The MHDM result is obtained from the 2HDM calculation |6^ by 
replacing |X| by XZ*: 

BR{B -^ Iv) = BR^^\B -^ £u)\l - i?|2, (58) 

where the SM result is 



BR^'^^B - iu) = -^'-^-^ I 1 - £1) mVubl'TB, (59) 

and 



Svr \ m^ 



In 2HDM (Model II) the interference term always reduces the rate, while in MHDM it may also enhance 
it, depending on the relative phase between X and Z. Nevertheless, in order to get bounds on \XZ\ 
one should take the minimal value for the branching ratio, which corresponds to a.rg{XZ*) = 0. Thus, 
the bound on \X\ in 2HDM is the same as the bound on \XZ\ in MHDM. 
We conclude that large contribution to B ^ iv may arise for large XZ. 
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Figure 6: The upper bound on \Y\ as a function of the lightest charged scalar mass mn- The three 
curves correspond to rrit = 100 (solid), 140 (dashed) and 180 (dotted) GeV. 

4 Combined analysis 

In this chapter we combine the results from the previous chapter and determine the region in parameter 
space allowed by the data. 

Several processes are sensitive to |y|. For some of them there are very poor experimental data: 
Bg — Bs mixing and B — > Xsi'i^- For others, the experimental bounds are far above the SM predictions: 
K"^ -^ vr+z^P, B — > i^i~ and B — > Xgi^i~ . The bounds obtained from Amx and Ki — > /x"*"/!" are 
very weak. The decay b ^ sj gives a strong constraint on |y| in the case of 2HDM {Model II) ||69|| . 
In MHDM, however, there is a possibility of cancellation between the two terms that depend on the 
MHDM parameters (see eq. (p6|)), and such a bound cannot be obtained. Then, there are three 
observables that constrain \Y\: e (eq. (|2l|)), Xd (eq. (|l7|)) and e;, (eq. (|3^)). At present, ef, gives the 



strongest bound |36]. It is important to note that the bound from e^ is strong since the experimental 
value is far from the SM prediction. Actually, at the la level mt is found to be less than the CDF 
lower bound. Had the data fitted large mt, the bound on |y| from B — B mixing may become the 
strongest. Our results are given in fig. g. For example, we obtain 

\Y\ < 1.3 for mi=140 GeV and mn = 45 GeV. (61) 



The only potential constraint on \X\ comes from D — D mixing (eq. (|29|)). However, the current 
experimental bound with the large uncertainties in the CKM elements gives no useful bound. We 
conclude that |^| is constrained by the perturbativity bound (eq. (|lq)): 

\X\ < 130. (62) 

\Z\ is constrained by e~fi universality in r decays (eq. (^)) and by the perturbativity bound 
(eq. (|l^)). For mn > 175 GeV the perturbativity bound is stronger: 

\Z\ < min(1.93 tuh GeV-\340). (63) 



\XZ\ is constrained hy B ^ Xtu-j- (eq. (^)), B -^ fj,u, B ^ ti/ (eq. (pq)) and by the product of 
the bounds on \X\ and \Z\. For light charged scalar B -^ Xtv^ give the strongest bound, while for 
niH > 370 GeV the perturbativity bound is the strongest: 

\XZ\ < min(0.32 mjj GeV"^, 44200). (64) 



In 2HDM {Model II), where X = Z, eq. (p^ ) gives \X\ < 0.56 rriH GeV , which is weaker than 
the published bound |^, 54 1. The difference arises due to the different values calculated for the 



interference term and the different values taken for the SM branching ratio. 
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Figure 7: The upper bound on \XY\ as a function of the lightest charged scalar mass mn- The three 
curves correspond to mt = 100 (solid), 140 (dashed) and 180 (dotted) GeV. 



Figure 8: The upper bound on Ini(Xy*) as a function of the lightest charged scalar mass mn- The 
three curves correspond to mt = 100 (solid), 140 (dashed) and 180 (dotted) GeV. 

It is important to note that below certain values of \XZ\ and mn the destructive interference leads 
to suppression and a bound cannot be obtained [68|. This happens at \XZ\ < 0.18 mj^ GeV^ for 



B -^ XtVt and at \XZ\ < 0.07 m^ GeV~ for B — > £u. We see that although the current bound from 
B -^ Xtv-t is stronger, B ^ iu can potentially give a stronger bound. 

The bound on Re(XZ*) is the same as the bound on \XZ\. However, the bound on the CP 
violating parameter lui{XZ*) is stronger. Im(XZ*) can be bounded from the same processes that 
bound \XZ\. For light charged scalar B — > XtUt give the strongest bound, while for mn > 440 GeV 
the perturbativity bound is the strongest: 

lm{XZ*) < min(0.23 mj^ GeV'^ 44200). (65) 

Several processes are sensitive to \YZ\. For some of them there are very poor experimental data: 
B -^ Xsfv, B -^ T^T^ and B -^ XgT^r^ . The only measurement that could potentially bound \YZ\ 
is K^ -^ ■K^uv (eq. (p9D). However, due to the current experimental bound and the large uncertainties 
in the CKM elements this bound is weaker than the product of the bounds on \Y\ and \Z\. The bounds 
on lva.{YZ*) and on Re(Xy*) are the same as the bound on \YZ\. For example, we obtain 

\YZ\ < 110 for mi=140 GeV and mu = 45 GeV. (66) 

The processes that are sensitive to XY* are B — B mixing, B -^ Xsi'^i~ and b -^ s^. The 
strongest bound comes from b ^ S'j (eq. (|3^)). The bound on \XY\ is very sensitive to arg(Xy*) and 
|y|. The weakest upper bound on \XY\ is obtained for maximal \Y\ and arg{XY*) = vr. Then, the 
cancellation between the two terms, the one proportional to XY* and the one proportional to |y| , is 
maximal. Our results are given in fig. ^. For example, we obtain 

\XY\ < 4 for mt=UO GeV and tuh = 45 GeV. (67) 

The bound on Jie^XY*) is the same as that on |Xy|. However, the bound on the CP violating 
parameter lm{XY*) is stronger. lni(XY*) is bounded by CP violating processes: NEDM (eq. ([50|)) 
and e (eq. ^4|)), and by CP conserving ones: B — B mixing, B -^ Xsl'^l~ and b -^ s^. The strongest 
bound comes from b ^ sj Q. Our results are given in fig. |8|. For example, we obtain 

Im{XY*) < 2 for mt=UO GeV and mn = 45 GeV. (68) 
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5 Discussion 

In this chapter we discuss the differences between the general MHDM and the widely discussed 2HDM 
{Model //)0. In 2HDM there are only two parameters in the charged Higgs sector: the mass of the 
charged scalar, niH, and the ratio of the vacuum expectation values, tan/5. The lower bound on tan/3 
is the same as the bound on 1/\Y\ obtained from the decay Z — > bb. Since in 2IIDM X = Z, the upper 
bound on tan [3 is the same as the bound on ^\XZ\. For intermediate values of tan /? the bound on 
ran is obtained from the decay b ^ sj (see e.g. fi^]). There are two main reasons for the differences 
between MHDM and 2HDM: 

• In MHDM there are four parameters in the charged Higgs sector. In 2HDM, XY* = 1 and 
X = Z so there are only two. Thus, processes that involve combinations of these parameters 
can be enhanced or suppressed in MHDM compared to 2HDM. 

• The bounds on the parameters in MHDM may be weaker than the corresponding bounds in 
2HDM. Thus, the possible effects are larger. 

5.1 Bound on mn 

In 2HDM b ^ sj gives a lower bound on mn almost independent of tan/3: rriH^WO GeV. In MHDM, 
however, this bound does not hold. This bound has implications for the experimental search of the 
charged scalar. In LEP200, a charged scalar of mass up to approximately 80 GeV can be found |19|. 



Consequently, the charged scalar of 2HDM cannot be detected at LEP200 while the charged scalar 
of MHDM can. Furthermore, the charged scalar of 2HDM decays dominantly to tv and cs. Within 
MHDM, when \X\ is large and \Y\ and \Z\ are small, the cb channel becomes important. Large 
background from W decays is expected if rriH ~ myy. Since the W^-boson hardly decays into cb, the 
H~^ — > cb decay mode is easier to detect. 

5.2 Processes that depend on |X| or \Z\ 



In 2HDM the bound on X = Z is the same as the bound on \J\XZ\ in MHDM. Consequently, the 
bounds on |X| and \Z\ are weaker within MHDM. Thus, processes that depend on \X\ or \Z\ may be 
enhanced within MHDM compared to the 2HDM case: 

• The charged scalar box diagrams can saturate the experimental bound on D — D mixing. Fur- 
thermore, measurable CP violating effects in the interference of mixing and decay may arise. 
Intuitively, the reason is that the charged scalar box diagrams introduce dependence on third 
generation parameters. Based on the discussion in p^ we find that the measured asymmetry 



_ N[{rX)n{P+p-)D]-N[{£+X)D{P+p-)n] 
^^ - N[{i-X)DiP+P-)D] + N[{i+X)D{P+P-)D] 



'in this chapter we refer only to Model II of 2HDM. 
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[P'^P is a CP eigenstate, e.g. tt+tt or K^K ) can reach its current upper bound 

Acp < 0.16. (70) 

Note, that no new source of CP violation beyond the phase of the CKM matrix is required. A 
recent estimate [|70|| shows that within one year operation of a Tau-Charm Factory asymmetries 
of order 1% are measurable. 

Bg — Bg mixing can be enhanced. The ratio between the two mixing observables, x^ and Xs, is 
given by: 



|2 

xd \Vtd\ 






Deviations from r = 1 in eq. (^) are due to flavor SU{?>) breaking effects. Prom eq. ( [l8|) we 
see that when |X| is close to its upper bound and rrifj to its lower bound, the last term in I^h 
dominates Xs but it is still small for Xd- Then, r is enhanced. We find that even r ~ 7 is possible. 

• Violation of e-/i universality in r decays is proportional to \Z\'^. Since the bound on \Z\ is weaker 
in MHDM compared with 2HDM, the violation can be stronger. 

5.3 Processes that depend on Y Z* or XY* 

While in 2HDM XY* = Y Z* = 1, in MHDM XY* and YZ* are, in principle, arbitrary complex 
numbers. Thus, processes that involve these combinations can be enhanced or suppressed in MHDM 
compared with 2HDM: 

• FCNC decays with third generation leptons in the final state. In such processes the charged 
scalar box diagrams may be important. They depends on Y Z* and on the lepton mass. Thus, 
in processes where leptons from the first and second generation are involved the box diagram is 
negligible. For K^ — > tt'^vv we find that the SM upper bound 

BR^*^(A"+ -^ -K+vu) < 4 X 10"^°, (72) 

is modified into 

which is more than 3 times weaker than the SM bound. For B -^ Xgvv we find that the SM 
upper bound 

BR^^\B -^ Xsvu) < 5.5 X 10-^ (74) 

is modified into 

which is more than a order of magnitude weaker than the SM bound. A similar enhancement 
may occur for the decays B — > t^t^ , Bg -^ t^t^ and B -^ XsT^t^ . 
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ruH GeV 


1^1 


\y\ 


z\ 


\XY\ 


\XZ\ 


YZ\ 


lm{XY*) 


lm{XZ*) 


45 


130 


1.3 


85 


4 


650 


110 


2 


465 


200 


130 


1.9 


340 


8 


12800 


650 


4 


9200 



Table 2: The upper bounds on the MHDM parameters for a representative value of rrit = 140 GeV. 

• Processes that involve the photon penguin diagrams depend on XY* . In 2HDM, the intermediate 
charged scalar enhances the effective bs'j vertex. In MHDM, the intermediate charged scalar 
may cause larger enhancement or, when aTg(XY*) ~ vr, suppress the effective bs'j vertex. Thus, 
processes like b ^>- s'j and B — > Xsi~^i~ may have different predictions in MHDM and in 2HDM. 

5.4 CP violation 

Within MHDM additional CP violating phases in charged scalar exchange are allowed. Those phases 
are absent in 2HDM. CP violation from charged scalar exchange can make only a very small contri- 
bution to e (less then 4%), and cannot be the only source of CP violation |Q. The charged scalar 
exchange may have a small effect on CP asymmetries in neutral B decays, at most 0.02 shift in the 
measured CP asymmetries Q. However, large contributions to the NEDM and to CP violation in top 
decays [^] are possible. 

6 Summary 

We have studied the charged Higgs sector of the general Multi-Higgs-Doublet Model (MHDM) . Using 
experimental data and the requirement of perturbativity we constrained the model parameters. These 
bounds are summarized in Table ^ for a representative value of mt = 140 GeV. 

We pointed out differences between the general MHDM to the widely discussed 2HDM (Model II): 

• The bound on the charged scalar mass is lower in MHDM. Thus, the charged scalar of MHDM 
are in the mass range accessible by LEP200. 

• The bounds on \X\ and \Z\ are weaker in MHDM. Then, charged scalar exchange can saturate 
the experimental bound on D — D mixing, Bs — Bg can be enhanced and strong violation from 
lepton universality is possible. 

• Large effects are possible in FCNC decays with third generation leptons. In these processes the 
charged scalar box diagrams are important only in MHDM. 

• There are differences in the predictions for processes that involve the photon penguin diagrams, 
namely, 6 ^ 57 and B -^ Xgi'^i" ■ 
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• In 2HDM there is no source of CP violation from the Higgs sector. In MHDM there are new 
CP violating sources but only when the charged scalars are not degenerate. Large CP violation 
effects, that arise from the CKM matrix, are possible in D decays. 

We conclude that MHDM may have large effects in various processes. Many of these effects do 
not appear in the 2IIDM. 

Acknowledgments 

I thank Yossi Nir for his help throughout this work, Miriam Leurer and Zoltan Ligeti for useful 
conversations. I am grateful to Donald Finnell for providing me a computer program that enabled the 
calculation of V^. 

A Proofs of relations 

In this appendix we give proofs of relations between the coupling constants of MHDM. 
Relation 1 

n 

J2XiY* = l (for^^/cD^). (76) 

i=2 

For ^ii ^ ^u, unitarity of U gives 

^ [/,,[/*, =0 ^ Y.u,dJL = -Udiui^ => ^'t; ': ^" = -1- (77) 

i=l i=2 UdlU^^ 

Using the definition of Xi and Yi (eq. (^)) we get eq. (76). The proof of the other two relations in 
eq. (13) is obtained in a similar way. 



Relation 2 



v"" 



5:ix,|2 = --1. (78) 



i=2 ^d 

Using eq. (^), eq. (g) and the unitarity of U we get 

The proof of the other two relations in eq. (S) is obtained in a similar way. 



Relation 3 Xi^Yj and Zj. {i,j,k E 2..n) cannot be simultaneously large (more than \/n — 1), unless 
there is at least one Higgs doublet that does not couple to fermions. 



21 



The proof is given for a general case where there are m types of massive fermions and n Higgs doublets. 
We denote the coupling of each fermion type by X'^, where k is the fermion type. We demand 

|Xf I > a. (80) 

Using eq. (^) we get 

n 2 

a'<\X^f<j:\Xh" = "--l => v'>ia' + l)vl (81) 

i=2 ^fc 



Summing over all fermions types we get 



mv 



> {a' + 1) ^ vl (82) 

k=l 



We have to distinguish between two cases. First, m = n and each fermion type couples exactly to one 
Higgs doublets. Then eq. (|82|) with eq. (||) lead to 

nv'^ > (a^ + l)v^ ^ a < ^/n-l. (83) 

On the other hand, if there is at least one Higgs doublet that does not couple to fermions, a is not 
bounded. 

B Formulae 

In this appendix we give the expressions of various functions that play a role in our study. 

B.l Neutral meson mixing 

All the integrals that we have to calculate are of the form 



for n = 0, 1, 2. We get 



^ g" ln(a) b" ln(6) 

{a - b){a - c){a - d) ^ {b - a){b - c){b - d) ^ ^ 

c" ln(c) d" ln{d) 



{c - a){c - b){c - d) {d - a){d - b){d - c)' 
The integrals in the text are 

hiyt) = J2[i,i,yt,yt) = - ^ + 7; 73^ (86) 

(1 - yt) (1 - yt) 
w ^ wi 1 ^ 2 , {l+yt)Hyt) 

hiyt) = Ji{i,i,yt,yt) = ^ + ^ ^3— ' 

(1 - yt) (1 - yt) 

h{xt,XH) = J2{l,XH,Xt,Xt) = 
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Xt 



+ 



xjjlnixH) 



h{xt,XH) 



h{yb,ys) 



{xt-XH)il-Xt) ' (l-XH)ixt-XHf 
Jl{l,XH,Xt,Xt) = 



+ 



xtjxt + xtXH - 2xh) In(xf) 
(1 - Xt) {xt - xh) 



+ 



XH ln(x//) 



+ 



{xl - xh) In(xt) 



{xt - Xh){1 - Xt) ' {l-XH){xt-XH)'^ {l-Xtf{xt-XH)^' 

'h{'i^A,ys,yb)- 



The SM integral can also be obtained from the general integral J. However, we keep the conventional 
notation [^1 



io{y) 



1 



4(1-2/)' 



1 + 



h{x,y) = ln{y/x) 



3y 



2y My) 
y My) 



4(1 -y) 



1 + 



i-y 



For the calculations of the matrix elements we use the Vacuum Insertion Approximation 



Vll 



(S|[dV(l -75)6][(i>(l -75)6]|S) = i/lms, 



Sll = {B\[dil - I5)b][d{l - j5)b]\B) 



niB 



6 \nrn, + m^ 



fB^B- 



(87) 



B.2 Z 



with 



B.3 6 



66 



S7 



Gh^(x) 
G/,(x) 



c 



a VL 



27rsin^6'vy('u| + ^i) 



^'L 



1 1-2/1 

- + -sm^6'iy, 



VR = -sin^ 6'iy. 



12(1 -x)' 



(7 - 5x - 8x2)(l - x) + 6x(2 - 3x) ln(x) 



6(1 -x)- 

1 - 8x + 8x^ - x^ 



[(3 - 5x)(l - x) + 2(2 - 3x) ln(x)] 



12x2 ln(x) 



(90) 



(91) 



B.4 K^ 



TV ' UU 



DsM{x,y) 



1 j xy /4: — y 



y-x\l-y 

4-x^2 



DzH{x,y) 

l7{xt,Xe,XH) 
h{xt,Xi,,XH) 



y — x\l — X 
xy_ l-y + My) 
4 (1-y)' 

J2(l,Xt,X<?,XH), 
Jl(l,Xt,X£,XH). 



ln(y) - 2x - ( 1 

3 

+ 1 + 



3x 



(1-x)^ 



I — yj 1 — X 
X ln(x) > , 



(92) 
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B.5 B -^ Xri^r 

The phase space factors for three body decay with only one massless anti-particle in the final state 
are |7|, || 

Fps{x,y) = 12/ {s-x'--y'^){l-s)^g{s,x^,y^), (93) 

Jix+yf S 

Ki-J/)' ds 



with 



ryi—y) rj^ r, 

F;«*(x,y) = 6/ -(s-x^fg{l,s,y% 

g{a, b, c) = y^a2 + b"^ + c"^ - 2{ab + ac + be). (94) 



We do not give the analytic expressions of these functions. Fpg can be found in ||72||, while we use a 
computer program to calculate F^^^. 

B.6 e— /i universality in r decays 



The phase space function for the interference term can be found in |57] 



ps 

B.7 B -^ £+£- 



Fif{x) = 1 + 9x - 9x2 - x^ + 6x(l + x) ^^^^^^^ ^-g^^ 



1 / 2>x Sx^ ln(x) 

- xH \ 

2\ 1-x (l-x) 



G'^'ix) = l[x+rr^ + ^^^:^^^], (96) 



Jim) = l^l±i^, 

2(1 -y)2 

T f \ 7 n T n^ ^ - yt - ytlnjyt) 

hivt) = J2(l,l,yt,0) = "2 , 

(1 - yt) 
T ( N xn n^ a;iln(xt) XHln(xH) 

h{i[Xt,XH) = J2{l,Xt,XH,0) = - + 



{I - Xt){xH - Xt) {1 - XH){xt - Xh) 
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